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Brain Imaging and Neuroscience 


I keep six honest serving men 
(They taught me all I knew); 

Their names are Vhat and Vhy and When 
And How and Where and Who. 


Rudyard Kipling 
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Who and When and How: MEG and EEG 



• Measure the electric potential (EEG) 
magnetic field (MEG) generated by neural 
currents 

• Who: Reasonable spatial resolution 

• When: Real-time measures of brain activity 

• How: Frequency-specific connectivity 
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Contents 


• Humble beginnings 

• From one to hundreds of channels 

• Insights from the Maxwell’s equations 

• Moving to source space 

• Are the new sensor technologies a breakthrough? 
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First MEG with EEG reference 



Cohen D, J. Appl. Phys, 1967; Science, 1968 
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First Real-Time Magnetoencephalogram 




Frequency (Hz) 



Matti Hamalainen 08/2018 


6 





























































































I believe the high-sensitivity detection system used 
here shows promise both as a clinical and as a 
research tool for studying the brain. 

Cohen, Science, 1972 
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Sam Williamson + and Lloyd Kaufman (NYU): 

• First VEF and mapping of somatosensory fields 

• Interdisciplinary approach 

• Magnetic Source Imaging (MSI) 
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Dipolar field pattern: Focal sources 


Somatically Evoked Magnetic Fields of the Human Brain 

Abstract. The human brain is found to produce a magnetic field near the scalp 
which iflrtrv i n s ynchrony w ith periodic electrical stimulation applied to a finger, Use 
of a highly sensitive superconducting quantum interference device as a magnetic 
field detector reveals that the bruin's field is sharply localized over the primary pro¬ 
jection area of the sensory cortex contralateral to the digit being stimulated. The 
phase of the response at the stimulus frequency varies monotonically with the repeti¬ 
tion rate and at intermediate frequencies yields a latency of approximately 70 milli¬ 
seconds for cortical response. 


We report here the detection of mag¬ 
netic fields associated with the How of 
electric current in the brain in response 
to electrical stimulation of the fingers* 
Weak magnetic fields resulting from vi¬ 
sual stimulation have previously been 
detected outside the scalp {1-4). In con¬ 
trast to the diffuse nature of the visually 
evoked potential (VEPl, which is con¬ 


ventionally measured with scalp elec¬ 
trodes, the visually evoked field (VEF) is 
located over the visual cortex (2). This is 
to be expected as the electric currents 
giving rise to the VEF flow w ithin the vi¬ 
sual projection areas of the brain, while 
accompanying weaker currents in the 
skin remain undetected. The VEF has 
proved to be a significant indicator of 
brain function since its latency is directly 
correlated with the reaction time of a 
subject when the spatial frequency of a 
stimulus is varied (J, 6). The confined lo¬ 
cation of the VEF and its correlation 
with a measure of human performance 
prompted us to search for neuromagnetic 
responses evoked by other stimuli. The 
neuromagnetic response evoked by elec¬ 
trical stimulation of individual fingers of 
the hand—the somatically evoked field 
(SEF)—is similarly found to be well lo¬ 
calized over the primary projection area, 
in this case on the contralateral hemi¬ 
sphere in the region SI for the stimulated 
digit. The observed features of the SEF 
indicate that the neuromagnetic tech¬ 
niques can provide unique advantages 
for medical research. 

A split-ring electrode, fashioned from 
a rubber grommet and two stainless steel 


Fig. I. (a) Neuromagnetic field patterns on the 
left hemisphere for an electrical stimulus at B 
hertz applied to the little finger of the nght 
hand. Contours of equal magnetic flux in¬ 
dicate the relative amplitude of response for 
0.9, 0.7. and 0.5 of the maximum response at 
the stimulus frequency, tbit The same pattern 
drawn on the conventional 10-20 electrode 
map Abbreviation: F.. fissure. 


Dipolar field patterns 
observed in several primary 
MEG responses: a current 
dipole seemed to be a 
reasonable model 


D. Brenner, J. Lipton, L. Kaufman, S.J. Williamson, 
Science, 1977 
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Lloyd Kaufman and David Cohen, October 2017 
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Soon SQUIDs Arrived to Finland 



Pekka Karp and Toivo Katila, 
carrying equipment from the wooden 
measurement cabin, around 1973 
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A Vision Can be a Game Changer 

whole-head MEG and 
10 fundamental reason 
uld not be built! 

Olli V. Lounasmaa 
Low Temperature Laboratory 
Helsinki University of Technology 

Finland 

We will establish a 
multidisciplinary team with 
full-time neuroscientists 

Riitta Hari, hired 1982 




MEG is not a hobby! 



We need 
there is i 
why it cc 
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An Early 7-Channel MEG Device 

Instrument Analysis Open-Source Software 


• 7-channel magnetometer 

• rf-SQUIDs, flux guides 

• 30-mm distance between channels 


Software on an HP 9825 
desk calculator 



0 : dsp "This 
p r o •? r o. m c o. 1 c u 1 a 
tes the lead 
field and Plots 
it f o r t. h e 
.seven " 

1 ! dsp "channel 
Magnet oc-et e r ? 
w h i c h m e o. s u r e s 
t, h e v e r t i c o. 1 
c o m p o n e n t o f 
2 : d s p " t h e m a g ri 
e t i c f i e 1 d c. t 
t h e ventices 
and at the cent 


3 : dsp “hexagon? 
t h e d i s t a n c e 
b e t w e e n t. h e 
n e a r e s t n e i s h b o 
u r s b e i n g 3 6 
n m " 

4 : d s p "th e Plot 


e t o m e t e r 
. The tr 


6 : dsp "to the 
different magne 
t o M e t e r s j n u m b e 
r 1 refers to 
t h e m a g n e t o n e t e 

7 : d s p " o' n t h e 
x-axis? n u m b e r 
2 66 degrees 
c o u n t e r c 1 © c k w i s 


.8: dsp ”Humber 
7 r e f e r s t o 
c en t e r. PIot 
i s n o r m a i i 2 e d 
so that m axi mum 

9: d s p " 1 e n g t h 
of 1 ine is unit 
y “ 

10: dim 917]? 

B [73 j C C 73 »Xt- 
1:73)YC-1:73 » 
ft* [13 

11: sol ~o ? 8>- 
1 2 ? 1 £ 

! £: 1.5 e - £ * D 5 
2EKXC13 j D*X [£3 g 
X C63 ; 8gY [ 1 3 g Y [4 
3-*Xt?3*V[?3 j 
f 3 D gY [23 ->V [33 

13: -0gX [33 gX [53 
;-2D*x [43;-rsog 
Y[53 gy[63 

14: 

-15: e n p “ t r a n s f o 
rmation vector" 
»fl[l] ?R[£3 j ri [33 
» fl [43 ? fi [53 f fi [63 


17: ggf-1 

18 : for 1=0 to 4 
19 : for J =0 to 4 
26 : 0 gS'gk 
21 : for K =1 to 7 
22: 1e-7((X [K3 — 

I / 10 8 ) 12 + (Y [ K 3 - 
J /10 0 ) T 2 + (Z x 
100 3 t 2 31 (- 1 . 5 )g 
F 

23 : F(Y[K 3 -J/ 
100 )gB[K 3 5 -F(X[ 
K 3 - 1 / 100 ) gCCK 3 
24 : fl[K 3 BlK 3 + 

SgS j R CK 3 C [K 3 + 
RgR 


r (R R + s S) g T 
if T>M? TgM 


:l: ent ’’ Want 
m a p ?" ? fi i 
!£•: if fli = "v" 


34 : ent X [-13 f 
X [03 ? Y [- 13 ? V ce 


t o X [ 0 3 

37: f o r J - V [-13 
to Y[ 83 
38: 8•*SgR 
39: 

40: for K=i to 7 
41: le-7it;X[K3- 
1/100)t £ +(YCK3- 
J/100)T2+(2/ 
180)t2)t(-1.5)g 
F 

42: F(V[K j-j• 
100 )gfi[K 35 -F(X[ 
K3-1/180)gCCK) 
43: fi tK3BCK3 + 

SgS j fi EK-3 C CK3 A 
RgR 

4 4: nex t K 
45: if RR+SS=0J 
Jmp 4 

46: Pit I-S/2M? 

J-R/2Mj1 
47: Pit I+ S/2 H ? 

J+R/2M> 2 
48: csi 2 .5 ?15 
1 b 1 

49: next J 
50: next I 
51: g t o 31 
52: end 
g11637 
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F.lmliolm GJ, Ilmoniemi RJ, and Wiik TO, Physica B+C, 107: 29-30, 1981. 13 






















































Neuromag 


Development of Multichannel MEG 

Systems in Helsinki 


SITRA 


Neuromag Oy 
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The recordings of magnetic fields over several 
brain areas simultaneously reveal qualitatively new 
information on the human brain and open up vast 
possibilities to investigate the interaction between 
separate brain structures in evoked studies and in 
the research of brain’s spontaneous activity. 

Ahonen et al., Physica Scripta, 1993 
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Dynamics of Cortical Activity 
in a Picture Naming Task 




■ sensory projection cortices 


nAm 
0 500 ms 


4 


— passive viewing 

— naming quietly 
naming aloud 


picture speech 

Hr 


Salmelin et al, Nature, 1994 
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Visual Stability During Eye Blinks 



Hari et al., Nature, 1994 
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Epileptic Spikes 




Forss et al., NeuroReport, 1993 
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“Statistics are really not necessary. I only conduct 

experiments in which the result is clear.” 

Ragnar Granit 
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The Sphere Model 
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Primary currents in the cortex 


No magnetic field from radial 
currents in the sphere model 







MEG = 0 


MEG * 0 
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current sources 








Auditory Responses 



“The distributions of the magnetic fields indicate that the auditory transient 
evoked response at a latency 100 ms as well as the auditory sustained 
response are generated at and around the primary auditory cortex.” 

Hari et a., Exp. Brain Res., 1980 
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Auditory N100 in EEG and MEG 


EEG 


t = 97.3 ms 
EEG step = 


EEG 



EEG step = 1.0 uV 



1.0 uV 



MEG 


t = 97.3 ms 
MEG step = 50.0 fT 



MEG 


t = 97.3 ms 
MEG step = 50.0 fT 
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Sphere model: Early contributions 


• Surface integral formulation for a piecewise 
homogeneous conductor 

B = j- [7p xV'-^rfV'+^Aff, f^V-IxA 

47r^ a k j K 

V 

Geselowitz, IEEE Trans. Magn., 1970 

• No magnetic field from radial dipoles 

• Layered sphere = homogeneous sphere 

• No contribution from volume currents to B, 

• Closed form expressions for the magnetic field 
using (vector) spherical harmonics as an initial step 

Arthur and Geselowitz, Biophys J., 1973 
Cuffin and Cohen, IEEE Trans. Biomed. Eng., 1977 
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We analyze the forward problem in the case of the 
spherical model, introducing simple ways of 
dealing with the tangential components of the 
magnetic field. 

Ilmoniemi, Hamalainen, and Knuutila, 1984 



Jukka Knuutila 
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Dipole in a sphere: Add radial sources 



Ilmoniemi, Hamalainen and Knuutila, 1984 
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Change to a homogeneous sphere 
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Summary of the sphere model graph 

theory 



Conductivity profile is irrelevant 



Simplified forward calculation 
A dry phantom can be constructed 
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The Homogeneous Head Model 



Skull and scalp taken into 
account 



Homogeneous model: 
skull taken as an insulator 


It is numerically shown that for the computation of B, produced by cerebral 
current sources, it is sufficient to consider a brain-shaped homogeneous 
conductor only. Comparisons with the sphere model are also included to 
pinpoint areas where the homogeneous conductor model provides essential 
improvements in the calculation of the magnetic field outside the head. 

Hamalainen and Sarvas, 1989 
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First MNE 


INTERPRETING MEASURED MAGNETIC FIELDS OF THE BRAIN: 
ESTIMATES OF CURRENT DISTRIBUTIONS 

Matti S. Hamalainen and Risto J. Ilmoniemi 

Low Temperature Laboratory 
Helsinki University of Technology 
SF-02150 Espoo 15, Finland 


REPORT TKK-F-A559 (1984) 


December 1984 

ISBN 951-753-362-4 

ISSN 0355-7790 

TKK OFFSET 




It would be desirable to improve the minimum norm estimates by finding ways to 
inject some a priori knowledge or assumptions of the experimenter. For example, 
one could confine the integration area to be the cortex, or one could require that 
the event in the brain be local. 

Hamalainen and Ilmoniemi, 1984 
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Retinotopic mapping with MNE 



Ahlfors et al. 1992 
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Are we still having fun? 


Jack Belliveau (MGH) 

• It is necessary and possible to combine multiple imaging 
modalities to benefit from their complementary strengths 
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Anders Dale 

• Computational models of anatomy based on MRI 

• First anatomically-constrained minimum-norm solutions 
( 1993 ) 
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Inflated Cortex 



The recovery of the topology of the cortical manifold makes it 
possible to conveniently view the solutions (as well as data from other 
activity imaging techniques) after gently “inflating” the brain, as well 
as allowing more effective intersubject comparisons. 

Dale and Sereno 1993 
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Inflation to a Sphere and Registration 


Individual 


Aligned with average brain 


Align sulcal patterns 

-► 

to the average brain 


MEG activity estimate 


Mapped to the average brain 


Morph 

-► 
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Source (brain) space instead of sensor space 



Sensor space map 



Source estimates 
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Visual percepts of an ambiguous scene 


MEG signals at an occipital sensor 



Percept ‘vase’ Percept ‘faces’ 
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Parkkonen et al., PNAS, 2008 
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Group analysis 


Significant activity in either tag frequency 



Amplitude ratio 
at ROI 


2.0 


CT> 

05 

+-> 

o d) 
o 

05 05 

5— M— 

0 ^ 

"O O) h c 
D (is 1.0 


Q. 0 

5 A 
< > 


1.0 


'faces' 'vase' 

Percept 



Parkkonen al., PNAS, 2008 





Left 


75% 


50% 

N = 
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Interacting parallel pathways associate 
sounds with visual identity in auditory cortices 


Object tracking task trials 


Visual identity and voice unchanged (AV Repeat ) 


+ 




Visual identity changes, voice unchanged (AV Change ) 



Effects of visual identity changes 
in auditory cortex 


Left superior 

temporal cortex 250-350 ms 



HG 


PT 


PP 


STG 


Cluster p 

^^Change > 

AV Repeat 

0.05 0.001 



Ahveninen J, Huang S, Ahlfors SP, Hamalainen M, Rossi S, Sams M, Jaaskelainen IP. Neuroimage, 2016. 124: p. 858-868. 
PMC4651767 







Different Frequency Bands in Connectivity 


y (50-90 Hz) 



* f 

cortex W f 4 




a (10 Hz) ll 

0 (7 Hz) \\ 

subcortical 

structures 

J (3 (20 Hz) 

Robert Oostenveld 


Matti Hamalainen 08/2018 


40 







Spatially and Temporally Specific Differences 

in Long-Range Connectivity 



• FFA seed, normalized (Z) coherence 

• Alpha band (8-12 Hz) 

Khan et al., PNAS, 2013 
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• Houses > emotional faces 
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Local Connectivity Different in FFA only 



• Phase-Amplitude Coupling differences between emotional faces and 
houses only in controls 

• Alpha band (8-12 Hz) phase and gamma band amplitude coupled 


Khan et al., PNAS, 2013 
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Modulation Index 




















































The Signatures Correlate wiith Behavior 


+ TD 
ASD 



a> 
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^-Coherence Precuneus 
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Small world property Small world property 


Cortical Beta and Gamma Rhythm Resting-State 
Networks Follow Distinct Maturation Trajectories 



10.5 14 17.5 21 24.5 28 
Age 


Khan et al., Neuroimage, 2018 
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Brainstorm is a collaborative open-source application dedicated to magnetoencephalography (MEG) and electroencephalography 
(EEG) data visualization and processing, with an emphasis on cortical source estimation techniques and their integration with 
anatomical magnetic resonance imaging (MRI) data. The primary objective of the software is to connect MEG/EEG neuroscience 
investigators with both the best-established and cutting-edge methods through a simple and intuitive graphical user interface 
(GUI). 
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Research Article 

FieldTrip: Open Source Software for Advanced Analysis of MEG, 
EEG, and Invasive Electrophysiologjcal Data 

Robert Oostenveld , 1 Pascal Fries , 1,2 Eric Maris , 1 and Jan-Mathijs Schoffelen 1 
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This paper describes Field Trip, an open source software package that we developed for the analysis of MEG, EEG, and other 
electrophysiological data. Ibe software is implemented as a MATLAB toolbox and includes a complete set of consistent and user- 
friendly high-level functions that allow experimental neuroscientists to analyze experimental data. It includes algorithms for simple 
and advanced analysis, such as time-frequency analysis using multitapers, source reconstruction using dipoles, distributed sources 
and beamformers, connectivity analysis, and nonparametric statistical permutation tests at the channel and source level. The 
implementation as toolbox allows the user to perform elaborate and structured analyses of large data sets using the MATLAB 
command line and batch scriDtine. Furthermore, users and develooers can easily extend the functionality and inmlement new 
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nAm 


Bridging 


MEG Source 
Estimates 



Example subject. Avg n=100 trials 

Steven Stufflebeam 
Matti Hamalainen 


Microscopic and Macroscopic 
Neurophysiology 


Computational Neural Animal 

Modeling Electrophysiology 



Predicted cellular level 
neural dynamics 

Stephanie Jones 


Cellular level details 
of neural dynamics 

Chris Moore 


Neural mechanisms of transient neocortical beta rhythms: 
Converging evidence from humans, computational modeling, 

monkeys, and mice, Sherman et al., PNAS, 2016 

MNE-Neuron Software: NIH R01EB02288 
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New Sensor Technologies 


• Improvements to conventional LowTc SQUID 
technology 

• HighTc SQUIDs 

- High sensitivity finally reached 

- Still needs cryogens 

• Optically Pumped Magnetometers (OPMs) 


Photodetector 


Vacuum Package 


Optical Fibers 



Monitor 

Photodetector 

Ferrule 


33 mm 
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Whole-Head MEG for Infants: the Baby MEG 



First whole-head infant MEG 
system with optimal sensor density 

Low-Tc SQUID system 

275 primary magnetometers 

35 three-axis compensation 
magnetometers 

Distance from sensors to room 
temperature (~ 6 mm) 

Cryocooler for closed-circuit He 
system 


• Ellen Grant & Yoshio Okada (BCH), Matti Hamalainen (MGH) 

• Supported by NSF 

Okada et al., Rev. Sci. Instrum., 2016 
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New Possibilities with New Sensors 


• On-scalp MEG 

• Adaptable MEG 

• Simultaneous MEG from 
several subjects? 

• MEG/TMS 
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Boto et al, Nature 2018 


ucyaiuiiem ot Microtechnology ( 




Conclusion: preliminary resu 
suggest that on-scalp MEG c 
distinguish phalanges on par 
better than in-helmet MEG 

Background: on-scalp MEG holds the 
promise for finer spatial discnmmair 

ssssssg 


Andersen et al, 2018 
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Comparison of Sensor Arrays 



500-channel On-scalp MEG VectorView 


■ 



500-channel EEG 
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Luessi et al., 2014 
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Localization Error (mm) Localization Error (mm) 


Localization Errors for Cortical Sources 





Luessi et al., 2014 
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Lifespan MEG: Premature babies to elderly 



• 64-channel flexible, conformal MEG (128 slots) 

• Optically pumped atomic magnetometers 

Svenja Knappe (PI): Conformal Pediatric Whole-Head MEG System with 
Optically-Pumped Magnetometers (5R01EB019440-03) 
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Integrated Multichannel TMS-MEG 


• 16 liquid-nitrogen cooled TMS probes 

• 25 atomic magnetometers for recording MEG 



Cryo-TMS probe 4 ms recovery of AM TMS-MEG array 

sensors after a TMS pulse 


Douglas Paulson (PI): Non-Invasive Integrated System for Brain Stimulation and 
Magnetocorticography (5 R44N S090894-03) 
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Summary 


• The multidisciplinary approach has been the key to 
advances in MEG instruments, methods, and 
applications 

• The advances in devices will inspire and necessitate 
advances in signal processing and data analysis 
methods 

• New breakthroughs are enabled by knowing and 
appreciating the past stumbling blocks and successes 
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Essentials 



Best people with different backgrounds work together daily: 

“Both information and gossip travel at the speed of light but 

information travels only 10 meters” 

“People are very stubborn there, very stubborn.” 

Molotov about Finland, 1974 
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Effect of sustained exposure to MEG 



Ole Jensen (20 y) 


Sylvain Baillet (23 y) Matti Hamalainen (37 y) 







1 > 1 


114 LJ 

' 
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“Sometimes statistics are really not necessary: 

three subjects is enough.” 
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